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Abstract 

A simple analytic theory for mixtures of hard spheres and 
larger polymers with excluded volume interactions is de- 
veloped. The mixture is shown to exhibit extensive immis- 
cibility. For large polymers with strong excluded volume 
interactions, the density of monomers at the critical point 
for demixing decreases as one over the square root of the 
length of the polymer, while the density of spheres tends to 
a constant. This is very different to the behaviour of mix- 
tures of hard spheres and ideal polymers, these mixtures 
although even less miscible than those with polymers with 
excluded volume interactions, have a much higher poly- 
mer density at the critical point of demixing. The theory 
applies to the complete range of mixtures of spheres with 
flexible polymers, from those with strong excluded volume 
interactions to ideal polymers. 

1 Introduction 

Athermal mixtures are remarkable because the only energy 
they possess is thermal energy; their behaviour including 
their phase behaviour is determined solely by entropy. But 
this does not mean that this phase behaviour is uninter- 
esting, entropy alone can drive phase transitions such as 
freezing and, in mixtures, demixing into coexisting fluid 
phases. This demixing is found almost whenever the two 
components of the mixture are very different, such as the 
athermal polymer and hard spheres considered here. We 
look at mixtures of hard spheres and athermal polymer 
where the polymer is as large or larger than the spheres: the 
mean-square end-to-end separation of the polymer greater 
than or equal to the diameter of the spheres. This com- 
plements earlier work on mixtures where the polymer is 
smaller. As might be expected we find extensive immisci- 
bility between the polymer and the spheres. We calculate 
phase boundaries and determine their scaling with the sizes 
of the components, for example, the density of monomers 
at the critical point is found to scale as one over the square 
root of the number of monomers in a polymer, just as does 
the critical point for demixing of a polymer and a poor 
solvent. 

We are most interested in the limit where the size of 
the polymer, Re, is much greater than that of the colloid, 



a. Re is the root-mean-square end-to-end distance of the 
polymer, and a is the hard-sphere diameter. Our theory 
is motivated by the idea that for fi^ > a the effect of 
the colloids is a small length-scale effect in the sense that 
if the colloid degrees of freedom of the mixture are inte- 
grated out we are left with polymers which at length-scales 
large in comparison to a behave qualitatively just like a 
polymer in structureless solvent. The colloids effect the 
quality of this effective structureless solvent but do not in- 
troduce any new physics on length-scales large with respect 
to a. If this is correct then the long length-scale behaviour 
and the phase behaviour, of colloid-polymer mixtures may 
be mapped onto the well-understood long length-scale and 
phase behaviour of a polymer plus solvent system. The 
phase separation of polymer and much smaller colloidal 
spheres will be qualitatively identical to that of a long 
polymer and a poor solvent. This phase separation has 
been extensively studied and is quite well understood [Q- 
[|] . The approximate theory we will develop is essentially a 
Flory-Huggins polymer-plus-solvent free energy with an ef- 
fective solvent quality which depends on the concentration 
of the colloid. It is known that a Flory-Huggins free energy 
describes rather well the qualitative features of the phase 
separation of a long polymer and a poor solvent; see the 
comparison with computer simulation data in Refs. |2f^|. 
It is also known that purely entropic effects can result in 
the polymer effectively being in a poor solvent Q . 

The literature on mixtures of colloidal particles and non- 
adsorbing polymers is extensive because colloid-polymer 
mixtures are common, and the limit we consider where all 
interactions are purely repulsive is a rather fundamental 
limit of these mixtures. Mixtures in which the polymer 
molecules are both larger than the particles and flexible (as 
opposed to semiflexible) are formed when the particles are 
small, a few nms across. Nanoparticles are colloidal parti- 
cles of this size as are proteins. Protein-polymer mixtures 
are common, for example polymers are mixed with proteins 
in order to induce the proteins to crystallise || ?j . Al- 
though in practice it is unlikely that the monomer-protein 
interaction is ever purely repulsive over all the surface of 
a protein, however the limit we study here provides a ba- 
sis for incorporating the effects of weak adsorption of the 
polymer onto a protein molecule. The limit in which all 
three interactions, sphere-sphere, monomer-monomer and 



Figure 1: A schematic of our mixture of large polymer 
molecules and colloidal spheres. The black discs represent 
the colloids and the curve represents a polymer coil. The 
rescaled monomers used to estimate the polymer-colloid 
interaction are drawn as dashed circles. 




sphere-monomer, are purely repulsive is an important and 
fundamental limit. It is fundamental in the sense that as 
we are assuming the monomers to be much smaller than 
the spheres, then the details of the monomer-sphere inter- 
action are irrelevant as long as it is purely repulsive; the 
details of the monomer-monomer interaction only effect the 
behaviour by altering a single parameter, the monomer- 
monomer second virial coefficient, and any sharply repul- 
sive sphere-sphere interaction will behave almost like hard 
spheres. Thus many of the details of the mixture are irrel- 
evant. 

The interactions between colloidal spheres and polymer 
molecules have been studied theoretically via a number of 
techniques: scaling approaches J^-fTo!, field theory Q-14], 
computer simulation |l6|| , integral equations |l7} 18], 
and other approaches |19|-|22|] . The phase behaviour has 
been studied via computer simulation |l5| |, scaling theory 
||, and perturbation theory p4[, |2^], w here Refs. 
are for ideal polymers. Reference |18| is a review of re- 
cent work on colloid-polymer mixtures, focusing mainly on 
the structure and on the results of integral equations, but 
also discussing other approaches and the phase behaviour. 
Earlier work by the author |^3| assumed that the phase sep- 
aration would occur when the polymer was semidilute. The 
more careful and better founded work here finds that this 
is not correct, the phase separation occurs at the boundary 
between the dilute and semidilute regimes where the free 
energy expression assumed in Ref. |23f| is not valid. There- 
fore, in particular the findings of Ref. |2^| for the polymer 
density where phase separation occurs are incorrect and 
should be discounted. The opposite limit to that of in- 
terest here, i.e., where the polymer molecules are smaller 
than the colloidal spheres, has been considered extensively, 
see Refs. ji~5|,|2(]-|29) and references therein. In this limit 
the polymer (with or without monomer- monomer excluded 
volume interactions) induces crystallisation of the spheres, 
there is no equilibrium separation into coexisting fluids. 



2 Model and physical picture 

Our colloidal particles are modeled by hard spheres of di- 
ameter er and our polymers are modeled by flexible chains 
of N monomers, each of length a. We characterise the in- 
teraction between a pair of monomers of the polymer with 
a second virial coefficient for this interaction, Bm- We 
start with a polymer in a good solvent far from the theta 
temperature where Bm — 0(a 3 ) and the polymer behaves 
as a self-avoiding walk (SAW). The interaction between a 
pair of hard spheres is similarly characterised by its second 
virial coefficient B$ = (2/3)ira 3 . The interaction between 
a sphere and a monomer is taken to be purely repulsive 
and short range, a monomer cannot overlap with a sphere 
but otherwise does not interact with it. The polymers do 
not adsorb onto the surface of the spheres. The monomer 
size a is much smaller than the diameter of the spheres, and 
here we are considering particles which are smaller than the 
polymers so our three length-scales satisfy a <C a < Re- 

A mixture of spheres and polymer is difficult to deal 
with because of the sphere-polymer interaction. Pure hard 
spheres and pure polymers have both been studied exten- 
sively and the equilibrium behaviour of both is well un- 
derstood. The sphere-polymer interaction is hard to deal 
with because of the large difference between the size of a 
sphere and that of a monomer and because of the con- 
nectivity of the polymer — if the interaction of a sphere 
and a monomer was independent of the interaction of the 
sphere with the next monomer along the polymer, then 
the sphere-polymer interaction could be estimated easily. 
However, when the monomers are much smaller than the 
sphere this is very far from being true. The problem of the 
disparity in sizes can be solved by rescaling the monomer 
size from a to a, i.e., by viewing the polymer not as being 
composed of N monomers of size a but of ub monomers 
of size a. This rescaling is quite common, see the book 
of de Gennes jjj. It has already been applied to mixtures 
of ideal large polymers and spheres p4[ , where it was re- 
ferred to as the extended Asakura-Oosawa model. We call 
the rescaled monomers of size a, blobs. In principle this 
rescaling can be done exactly, i.e., it can be done leaving 
the free energy etc. exact, but here we view it as part of 
a physically motivated approximation scheme. See Fig. [I] 
for a schematic of our sphere-polymer mixture showing the 
polymer composed of a chain of blobs of size a. If there 
are Nb monomers of size a in one blob of size a then ub 
is related to N by ub — N /Nb- So, relating the number 
of blobs to the number of monomers requires estimating 
Nb, we defer this to section [|. Until then we specify the 
polymer size by specifying n b and do not concern ourselves 
with how riB is found for a given polymer. The approxi- 
mate free energy we will obtain depends on the polymer 
length only through tib- 

Having performed the rescaling we approximate the in- 
teraction between a sphere and a polymer as being ns in- 
dependent sphere-blob interactions. This is quite a rea- 
sonable approximation as the blobs and spheres are of the 
same size. Applying this approximation before rescaling, 
i.e., approximating the sphere-polymer interaction by N 



sphere-monomer interaction is qualitatively wrong because 
the monomers are so much smaller than the sphere so 
the sphere interacts with many monomers at a time. We 
rescaled the monomers just so we could apply this simpli- 
fying approximation. The blob-blob and sphere-blob inter- 
actions are characterised by their second virial coefficients, 
Bb and B$b respectively. The second virial coefficient for 
the interaction between a sphere and a blob is of order that 
for the interaction between two hard spheres of diameter a 
but is a little smaller. We defer its estimation to section 
^. Returning to our assumption that the interaction be- 
tween a polymer and a sphere consists of ng independent 
blob-sphere interactions, this implies that the second virial 
coefficient for the sphere-polymer interactions is ubBsb- 
Finally, we remark that after rescaling, our mixture of poly- 
mers of blobs and spheres resembles the athermal polymer 
+ solvent mixture considered by Frenkel and Louis || . The 
phase separation in both cases is driven by unfavourable 
excluded volume interactions. 

3 Flory-Huggins— type theory 

We start with the basic Flory-Huggins theory for a poly- 
mer of monomers which interact only via excluded volume 
interactions, see for example the book of de Gennes, pages 
113-115. This is often referred to as an athermal polymer 
solution. For our 'monomers' we take the blobs not the 
original monomers of length a. Thus the 'monomer' density 
in the theory is actually the density of blobs which is equal 
to PpUb'-, Pp is the number density of polymer molecules. 
As usual we use a reduced 'monomer' density, </>, which we 
obtain by multiplying the blob number density by the vol- 
ume one blob excludes to another, 2Bb, 4> = ppub{2Bb)- 
The Flory-Huggins Helmholtz free energy F then has the 
usual form 

F(2B B ) 
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which defines the reduced Helmholtz free energy per unit 
volume /. Throughout, we use units such that the thermal 
energy kT = 1. This is for a polymer solution, no colloidal 
spheres present. We add on the contribution of the colloidal 
spheres using a virial expansion, 

/ = — In 0+ (1-0) ln(l -0) + 

(2B B ) {pc [hi Pc - 1] + Pc B s + PcRpnBB SB } , 

(2) 

where we have truncated the expansion after the second 
virial coefficient terms. We have dropped cubic and higher 
order terms, which is only valid at low colloid densities. 
The last term within the braces is the second virial co- 
efficient term for the polymer-sphere interaction: within 
our approximation it is just ns independent sphere-blob 
excluded volume interactions, each with an excluded vol- 
ume Bsb- In order to obtain a simple analytic theory we 
will neglect not only all terms for the sphere-polymer in- 
teraction beyond the leading order, second virial coefficient 



term, but all the terms from sphere-sphere interactions, in- 
cluding the leading order p c Bs term shown in Eq. (||). 
This latter approximation is quite severe but we do this 
in the expectation that the sphere-sphere interactions will 
not be very large when the colloidal spheres and polymer 
demix and that the sphere-sphere interactions, unlike the 
blob-blob interactions, are not essential to understanding 
the basic physics of this demixing. 

Making these two approximations and defining a reduced 
density of spheres (f>c — Pc{2Bb), Eq. (||) becomes 



/ = — ln0 + (1 - (j>) ln(l 
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where we have changed the In pc term to In <pc term which 
we can do as the difference between the two is a constant, 
\h(2Bb), which has no effect on the phase behaviour. The 
quantity b is the ratio between the excluded volumes of the 
sphere-blob and blob-blob interactions, b = Bsb/{2Bb)- 

Apart from a somewhat more complex dependence on cf>, 
Eq. (Q) is of the same form as the free energy of Eq. (1) 
of Ref. 0. Although in that work the mixture was of 
a mixture of thick and thin hard rods and the free en- 
ergy was exact. Below, we will transform Eq. (^) following 
the same approach as used in Ref. p0| . The free energy 
Eq. (||) is linear in the density of spheres, (pc- This makes 
it easy to analytically transform from fixed pp and pc to 
fixed pp and pc, where pc is the chemical potential of 
the colloidal spheres. This transform is useful as then we 
have a thermodynamic potential, called the semigrand po- 
tential, which depends on one density variable, pp or </>, 
and one field variable, pc- This is completely analogous 
to the Helmholtz free energy of a single component system 
in which temperature is important; that free energy also 
depends on a density variable (the number density) and a 
field variable (the temperature). As such once we have the 
semigrand potential calculating phase equilibria is just as 
easy as for a single component system. 

So, at fixed (f> and pc the relevant thermodynamic func- 
tion is the semigrand potential il. In fact it is slightly 
more convenient to work with the activity of the colloid 
zc — exp(pc) not the chemical potential. We need the 
semigrand potential f2 which is a Legendre transform of 
the Helmholtz free energy 
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which defines the reduced semigrand potential per unit vol- 
ume uj fn]] . The chemical potential pc is just the 4>c 
derivative of /, so, taking this derivative of Eq. (||), 



pc — In <f>c + 4>b = In Zc, 



(5) 



which can be rearranged to obtain an equation for <fic in 
terms of zc 



epe = zc exp (- 



(6) 



Using Eq. (||) to substitute for / and Eq. (||) to substitute 
for fi c , in Eq. (0), 



= -Mn0+(l-0)ln(l-0)-0 c , 
n B 



(7) 



but we want it in terms of the relevant variables which are 
4> and fic so we use Eq. (Q) to substitute for 4>c 

lj = — In <j> + (1 - <j>) ln(l -<j>)-z c cxp (-<j)b) . (8) 
n B 

This equation completely describes the thermodynamics of 
the mixture. 



3.1 The critical point 

We begin the determination of the phase behaviour by find- 
ing the critical point for demixing. First we need the chemi- 
cal potential \x of the polymer, which is just the (f> derivative 
ofEq. (fh, 
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+ zcbexp(-<j>b). (9) 



Now, the critical point is the point where the first and 
second <fi derivatives of the chemical potential are equal to 
zero. The derivatives are 



Oft 



1 



1 



n B (, 



1-0 
1 



n^ 2 (1 



- z c b 2 exp(-#) (10) 
4-zc6 3 exp(-<?(>6). (11) 



Setting them to zero results in two simultaneous equations 
for the polymer blob density c/) cp and the sphere activity 
z C q at the critical point. Combining these two equations 
yields an equation solely in terms of 4> cp , 

- (1 - (f> cp ) 2 + n B (<f P) 2 + - ^) 2 + 



bn B (^f (1 - ^P) 



0(12) 



and also one for 
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exp(0 c "6). (13) 



Equation (|l2| ) may be solved numerically for cjf p and then 
z°q obtained from Eq. (|T^). 

For large n B the equations simplify and we can solve the 
equations explicitly. For large n B we look for a solution 
with (j> c P small. Eqs. (fil) and (ph then yield 
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Figure 2: The variation of the volume fractions of polymer 
blobs and colloidal spheres at the critical demixing point, 
as a function of polymer length n B . The solid and long- 
dashed curves are for polymer in a good solvent, b = 3.9; 
the solid curve is the volume fraction of polymer blobs, 
rjgP and the long-dashed curve curve is that of the spheres, 
rf v . The dotted and dot-dashed curves are for polymer in 
a quite poor solvent, b = 20; the dotted curve is r\ c B and 
the dot-dashed curve is rf v . 
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In order to estimate the magnitude of the sphere-sphere 
interactions we require the volume fraction r\ of the spheres, 
V = /°c(7r/6)o' 3 = Pc(Bsb/8). At the critical point the 
volume fraction of the spheres is 
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n B > 1. (15) 



For large polymers the volume fraction r\ of spheres at the 
critical point is close to 1/86. For our neglect of the sphere- 
sphere interactions to be valid r\ must be small. We esti- 
mate 6 in the next section and find it to be generally around 
4 or larger, so the volume fraction of spheres at the criti- 
cal point is roughly 0.03 or less and our neglect of sphere- 
sphere interactions is not unreasonable. Having defined the 
volume fraction of spheres we can also define an effective 
volume fraction of the blobs. If we regard each blob as 
filling a spherical volume of diameter a then the 'volume 
fraction' of blobs equals 0(6/8) and so at the critical point 
we have a blob volume fraction at the critical point 



cp n B U „ i 

Tjg = §z n B > 1. 
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For large polymers, n B 3> 1, at the critical point the den- 

— 1/2 

sity of polymer blobs scales as n B while the sphere activ- 
ity tends to a constant as n B increases. The reduced den- 
sity of spheres (j>c tends to 6 -2 for n B large, from Eqs. ^) 
and (lil). 



^ In order to illustrate the trends in the demixing be- 
haviour with polymer length, measured by n B , we have 
calculated (using Eqs. ( |l2"|) and (|l^)) the volume fractions 
of spheres and of blobs at the critical demixing point, rf v 
and 77^, respectively. The results are shown in Fig. || where 
the solid curve is the volume fraction of blobs, and the long- 
dashed curve is the sphere volume fraction. For large n B , 
we see that while the sphere volume fraction is tending to- 
wards a plateau, that of the polymer blobs is continuing 



Figure 3: The phase diagrams of two colloid-polymer 
mixtures in the t\-t\b plane; the x and y axes are the vol- 
ume fractions of the colloidal particles and blobs, respec- 
tively. The curves denote the coexisting densities and the 
black circles denote the critical points. Both curves are for 
ub = 5 blobs. The solid curve is for a good solvent b = 3.9 
while the long-dashed curve is for a rather poor solvent, 
b = 20. The dotted lines are tie lines, lines connecting two 
coexisting phases. 
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interactions within a blob are strong, see the next section 
for a discussion of this point. The third virial coefficient is 
of order a 6 at the critical point, again assuming the intra- 
blob excluded volume interactions are strong. This is just 
as in the standard Flory-Huggins free energy for a polymer 
in a poor solvent Q. 

Thus, within our simple mean-field theory, adding small 
colloidal particles, a <C Re, to large polymer molecules in 
a good solvent is essentially equivalent to altering, wors- 
ening, the solvent quality. As more and more spheres are 
added the rescaled monomers, the blobs, start to attract 
each other and so the polymer and spheres demix just as 
polymer and a poor solvent demix. If this picture is correct 
then at length-scales large in comparison to the blob size 
a the polymers will behave just as a normal polymer in a 
good, theta or poor solvent, depending on the concentra- 
tion of spheres. The behaviour of polymers as the solvent 
quality is worsened and phase separation occurs has been 
well-studied and is now well understood; see the results of 
recent computer simulations These simulations have 

— 1/2 

found that the mean-field theory prediction of n B scal- 
ing of the density of polymer at the critical point and the 
size of the second virial coefficient Bi, are almost correct, 
there are only logarithmic corrections. 



to decrease, which is just what we expect from Eqs. ( |15[ ) 
and (|l^). For the calculations we have set 6 = 3.9 which is 
approximately its value in the a/a — > 0, Bm limit. We 
will discuss the estimation of b and its variation with sol- 
vent quality in section |4|. The phase diagram in the plane 
of the volume fractions of colloidal spheres and blobs, for 
spheres and polymers with ng = 5 blobs is plotted in Fig. ^|. 
Note the large region of fluid-fluid coexistence. 

3.2 Virial expansion 

We can expand out our semigrand potential Eq. (||) as a 
virial or density expansion 

ui = —\n<j) - z c - <j> + zbcj) +]- (l - z c b 2 ) <f> 2 + 
n B 2 

i(l + z c7 6 3 )^ + l(l-iz c 6 4 )0 4 + .... 

(17) 

So, the second B2 and third B3 virial coefficients for the 
monomer-monomer interactions in the presence of the col- 
loidal spheres are 



B 2 = \ (1 - z c b 2 ) (2B B ) 



z c b 3 ) {2B B f 



(18) 



Using Eq. ( |l4| ) for the activity of spheres at the critical 
point in the above equation for the virial coefficients, we 
find that the critical point occurs when the second virial 
coefficient in the expansions is negative and small, of order 
BbtIc, 1 ^ 2 ~ cr 3 ?^^ 1 ^ 2 , when is large. This holds for 



Bb of order er 3 which is correct when the excluded volume 



4 Solvent quality 

The phase behaviour depends on only two parameters: the 
number of blobs of size a in the polymer, n B , and the 
ratio of the blob-sphere to blob-blob excluded volume, b. 
The blob-blob excluded volume, 2Bb, is needed to convert 
from our reduced densities to number densities, but only 
the ratio b = B$b / (2Bb) affects the nature of the phase 
behaviour. 

First, let us consider polymers in which the excluded 
volume interactions are strong, we are far from the theta 
temperature, and the particles are not too small. This 
is the 'excellent' solvent regime of Odijk |2~l}| . Note that 
whether or not a solvent is excellent in this sense depends 
not only on the properties of the solvent but on the size 
of the particles. There is a parameter which describes how 
strong the excluded volume interactions are B,B2l, it is 
often denoted by C- For a single blob we have £b which is 

an 



a 6 
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(19) 



When (b -C 1 then an individual blob is close to being a 
random walk — the excluded volume interactions within 
a single blob are negligible. In the other limit, (b 3> 1 
the excluded volume interactions within a single blob are 
strong. 

The Cb — ► °o limit is the limit of a large SAW, as Nb — > 
00, <^b - > 00, Eq. (|H^). This is the scaling regime of an 
SAW, which is widely studied and employed, it is the limit 
in which Re scales as a power law of N with the Flory 
exponent v which is close to 3/5. Calculations on SAWs 
|32, 33 give the second virial coefficient between two SAWs 



Figure 4: The variation of the volume fractions of polymer 
blobs and colloidal spheres at the critical demixing point, as 
a function of solvent quality, measured by b. The solid and 
long-dashed curves are for a polymer of length jib = 3; the 
solid curve is the volume fraction of polymer blobs, rf^ and 
the long-dashed curve curve is that of the spheres, rf v . The 
dotted and dot-dashed curves are for a polymer of length 
riB — 30; the dotted curve is rfg and the dot-dashed curve 
is rj cp . 
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with radii of gyration Rq as Bb = 5-5Rq in the large 
(b limit. They also find R E = 2.5R G , so B B = 0.35cx 3 
between two SAWs with mean end-to-end separations of 
a. Hanke et al. [|l2j have applied field theory to obtain 
the result that BgB = 2.7 a 3 J34J. This result is not exact 
but is more than accurate enough for the purposes of the 
present theory. This result is obtained in the scaling limit of 
strong excluded volume interactions. For comparison, for 
ideal polymers B$b = 3. Oct 3 . A swollen blob is more open 
than one which is ideal so B$b is correspondingly smaller 
for polymers with the same Re- So, with strong excluded 
volume interactions even within blobs the ratio b = 3.9. In 
this limit the blob-blob interaction is strong and the chain 
as a whole will be swollen so if the radius Re of the chain is 
known then jib may be estimated from iib ~ (Re/ct) 5 ^ 3 . 
Alternatively, if the number of monomers Nb in a chain 
with R E — a is known then the number of blobs may be 
found from ns = N/Nb- 

So far we have considered only polymer-colloid mixtures 
in which the solvent for the polymer is sufficiently good 
and the colloidal particle sufficiently large that (b ^S> 1 
and even pieces of the polymer as small as a are strongly 
swollen. Then we can use the value of the blob-blob and 
blob-sphere second virial coefficients, Bb and B$b, in the 
Cb — > co limit. But what if the solvent quality is less good 
and the particles not too large? As the solvent quality de- 
creases the monomer- monomer interaction, Bm decreases 
from its value in a good solvent which is of order a 3 . This 
will decrease the blob-blob interaction, measured by Bb, 
while leaving the blob-sphere second virial coefficient B$b 
still at around a 3 . The second virial coefficient for the 
interaction between an ideal chain, Re = <r, and a hard 



sphere, of diameter a, is known exactly [ pT[ and is close 
to 3cr 3 . Thus, when £ B is no longer much larger than one, 
the ratio between the sphere-blob and blob-blob excluded 
volumes, b, will increase as C,b decreases. It is divergent 
for ideal polymers as then the blob-blob excluded volume 
is zero. 

We expect that as the solvent quality for the polymer 
worsens and the polymer-polymer interactions weaken that 
phase separation will be enhanced, the polymer and colloid 
will be less miscible. With this in mind we return to the 
equation for the density of polymer at the critical point, 
Eq. ([[2]). We assume that the critical density will be very 
low but make no further assumptions, Eq. (|l^) then sim- 
plifies to 



-1 + n B (0 cp ) 2 + b(f> c P + bn B (<j> cp Y 



0. 



which is a quadratic solution with a physical root 

, cp= -1 + Vl + 4(n 5 /b)(l + l/b) 
9 2n B (l + 1/6) 



(20) 



(21) 



For ris/b large, i.e., long polymers with blob-blob interac- 
tions which are not too weak, this equation simplifies to 
Eq. (^) — for very long polymers even weak monomer- 
monomer and hence blob-blob excluded volume interac- 
tions are sufficient to obtain the standard scaling of the 
critical density as n^^ 2 . The other limit in Eq. ( pl| ) is when 
ns/b is small, i.e., the blob-blob interactions are weak and 



the polymer is not too long. In this limit, Eq. (21 ) simplifies 
to 



n B /b<. 1, 



r& p = l/8 n B /6«l. 



(22) 



(23) 



The density of blobs at the critical point is independent of 
the length of the polymer. Recalling that 4> is the number 
density of blobs times 2Bb, we rewrite Eq. (p2h as 



pp = 



1 



n B B, 



SB 



n B /6< 1, 



(24) 



where pp is the number density of polymer molecules at the 
critical point, which is of order l/(ns(T 3 ). This scaling has 
been derived before for mixtures of colloidal hard spheres 
and ideal polymers [24f| . Our present theory for polymers 
with excluded volume interactions, although rather differ- 
ent in a number of ways from the theory of Ref . |2ij] predicts 
the same critical density of polymer when these excluded 
volume interactions are turned off. Giving us confidence in 
both theories. We can easily obtain the density of spheres 
at the critical point in the ns/b <C 1 limit. The reduced 



activity is z C q = e/(ns6), from Eq. (y 
reduced density at the critical point 
finally the volume fraction 



b c S ^ 



and hence the 
1/(71b&) and 



1 



8n B 



n B /6< 1, 



(25) 



which again is the same scaling as found previously p4j . 
For weak interactions between the blobs the density of 
spheres at the demixing critical point is dramatically be- 
low its value when there are strong interactions between 
the blobs, Eq. (|l5|). Also, note that the ( parameter for 

the whole chain is ( P = N 1 / 2 B M /a 3 ~ n^ 2 B B /cT 3 . We 

1 /2 

can rewrite this in terms of b, £p ~ n B /b, ignoring a nu- 
merical prefactor. When ris/b <C 1, then substituting for 

— 1/2 

b, £p <C n B and so as is larger than one £p 1 — an 
individual polymer molecule in dilute solution and before 
the colloidal spheres are added is ideal, its Re — aN 1 / 2 . 
Equations (22) to (p5|) are for demixing of a polymer which 



is effectively ideal, its monomer- monomer excluded-volume 
interactions are negligible. 

Returning to Fig. ^, we can compare the densities at the 
demixing critical point for polymers in a good solvent (solid 
and long-dashed curves) with those in a rather poor solvent, 
b = 20 (dotted and dot-dashed curves). We see that as 
the solvent quality worsens the volume fraction of polymer 
blobs increases (compare the solid and dotted curves) and 
the volume fraction of spheres decreases (compare the long- 
dashed and dot-dashed curves). Also, for b — 20 and for 
not-too-large ng the density of blobs at the critical point for 
demixing is relatively insensitive to Tip, which is what we 
expect from Eq. ( p3| ) . Note that for the polymer in a good 
solvent the effective volume fraction of polymer blobs is 
lower than that of the spheres at the critical point whereas 
for a poor solvent, b = 20, the opposite is true. 

In Fig. ^ we have plotted the the variation of the den- 
sities at the critical point for fixed ns and varying sol- 
vent quality b. For a small polymer, n E — 3, the vol- 
ume fractions of blobs and spheres are plotted as the solid 
and long-dashed curves, respectively, while the dotted and 
dot-dashed curves are the volume fractions of blobs and 
spheres, respectively, for a much longer polymer, ub = 30. 
The mixture of the longer polymer and the spheres demixes 
at lower densities than the mixture of the shorter polymer 
and spheres of course, and as the solvent quality decreases 
the volume fraction of polymer increases while that of the 
spheres decreases. In Fig. y we have plotted the phase dia- 
grams in the plane of the two volume fractions for polymers 
of length ns — 5 for a good solvent, b — 3.9 and a rather 
poor solvent, b = 20. The coexistence curve for the poorer 
solvent lies outside that for the good solvent: reducing the 
solvent quality increases the extent of the immiscibility. 

The larger density of polymer at the critical point comes 
from the fact that, near their critical points, the third virial 
coefficient is relatively much larger for ideal polymers than 
for polymers with excluded volume interactions. The val- 
ues of the virial coefficients are obtained by inserting the 
activity of colloidal spheres at the critical point, 



cp 
Z C ' 



into 



Eq. (|l^). For polymers with excluded volume interactions 
near the critical point B3 — 0(a 6 ). While for ideal poly- 
mers near their critical point B3 = 0(cr 6 /ns)- Small B^'s 
lead to high critical densities, see the appendix for details. 



4.1 Comparison with computer simulation 

Finally in this section, we compare with the results of re- 
cent computer simulations of colloidal particles and hard 
spherical particles by Meijer et al. |55). They studied col- 
loidal particles + SAWs with R E /a = 4.8, 7.0 and 9.9. In 
the simulations the size of colloidal particle not that of the 
polymer was varied but for the moment we will assume that 
in each case ^> 1 so that in the simulations the blob-blob 
and blob-sphere interactions are both in the good solvent 
scaling regime and our parameter b — 3.9. This leaves us 
with the problem of estimating the values of our parameter 
Ub for the simulated systems. For large polymers we must 
have that tib ~ {Re/ct) 1 / 1 ' as ub is extensive in the con- 
tour length of the polymer. In order to obtain an estimate 
for ub we set the unknown numerical prefactor in this scal- 
ing relation to 1, and v = 0.6, and so obtain = 14, 26 
and 46 for the three simulated systems. 

For these three systems we predict critical points at poly- 
mer blob volume fractions rfg = 0.045, 0.035 and 0.028, 
respectively. Preliminary simulation results for the critical 
densities are p C p '/ 'pp = 2.04, 3.19 and 4.65, where pp = 
l/[(47r/3)i?^.]. Using the theoretical result R E = 2.5R G 
in the good solvent regime and converting from Re/ct to 
ns as above we have that pp — 3.7 /n B ^. Converting our 
theoretical predictions to values of the ratio p C p/pp, we 
have p C p I pp = 0.19, 0.25 and 0.31. The theoretical predic- 
tions are about an order of magnitude too small, although 
the trend with increasing Re/ct is correct. Part of the 
discrepancy may come from our crude estimation of the 
relationship between the size of the polymer, Re/ct, and 
the number of blobs, tib, but it seems very likely that the 
theory is also underpredicting the density of polymer when 
it demixes from the spheres. Preliminary results for the 
volume fraction of spheres at the critical demixing point, 
i] cp , are around 0.2-0.25 whereas we find 0.03-0.05 for this 
range of values of Re/ct. Again the theory overestimates 
the extent of the immiscibility. It should be noted that, 
in simulation, for the largest value of Re/ct, the colloidal 
spheres have a diameter only about ten times that of the 
monomer. If the blobs are too small to be in the good sol- 
vent regime, the requirement Cb 3> 1 will be violated, and 
our parameter b will be > 3.9. Then the weaker blob-blob 
interactions in simulation will increase the polymer density 
at demixing with respect to that given by the theoretical 
prediction for the good-solvent scaling regime. The sim- 
ple theory derived here is clearly not quantitative but this 
is perhaps no surprise, it is really only capable of giving 
rough estimates and the qualitative nature of trends. The 
precise nature of the trends for Re/ct ^> 1 will be those for 
a long polymer in a poor solvent; see the simulation results 
of Refs. §%. 



5 Conclusion 

Mixtures of hard spheres and larger, flexible polymers 
which do not absorb onto the surface of the spheres ex- 
hibit extensive immiscibility. The cross excluded volume 



interactions in these mixtures, i.e., the excluded volume in- 
teraction between the sphere and the polymer, are large. 
So, the spheres and polymers 'get in each others way' so 
reducing each others entropy and driving them apart into 
separate phases || . The tendency to demix increases as the 
polymers become larger and larger, because the excluded 
volume interactions scale linearly with the length of the 
polymer while the translational entropy gained by mixing 
solutions of spheres and polymers does not vary with poly- 
mer length. Here we derived a simple analytic theory for 
these mixtures and found that when the excluded volume 
interactions were so strong that even blobs of size equal to 
that of the spheres were swollen, a mixture of polymer and 
much smaller spheres behaves much as a polymer in a poor 
solvent does. By contrast, as we showed in earlier work, a 
mixture of spheres and ideal polymers behaves rather dif- 
ferently. There the effective third virial coefficient of the 
polymer is very small, which pushes up the polymer density 
at the critical point for demixing. 

Our results are of relevance to mixtures of globular pro- 
teins and polymers, as in these mixtures it is easily pos- 
sible to have polymers large than the protein. However, 
our assumption of a purely repulsive interaction between 
the polymer and the spheres is rather unrealistic for pro- 
teins which have rather complex surfaces. Some part of 
this complex surface may well attract the monomers. A 
clean comparison with experiment could be done with ex- 
periments on either small synthetic colloidal particles or 
nanoparticles. Another possibility is to instead of making 
the colloidal particles smaller, make the polymer bigger by 
employing DNA |56|. There however, the colloidal particles 
would have to be reasonably large as the effective monomer 
length of DNA is about a = lOOnm and our theory assumes 
that a a. An objective for future work could be to relax 
this restriction to account more accurately for mixtures of 
spheres with semiflexible polymers such as DNA. 

It is a pleasure to thank P. Bolhuis, A. Louis and E. 
J. Meijer for inspiring conversations, and for sending me 
results of their computer simulations, and A. Hanke for 
providing me with his calculated value for B$b- I would 
also like to thank D. Frenkel, also for inspiring conversa- 
tions, and for the invitation to visit AMOLF, where this 
work was started, 

Appendix: Virial expansions and 
critical densities 

Here we explore how the critical density varies with the size 
of the third virial coefficient. Consider the simplest possi- 
ble virial expansion that has a critical point, an expansion 
truncated after the third virial coefficient. For the chemical 
potential /i as a function of density p this is 

H = hip + 2B 2 p + 3B 3 p 2 , 

where B 2 and B 3 are the second and third virial coeffi- 
cients, respectively. Assuming that B 3 is fixed and that B 2 
varies linearly with some temperature-like variable t, we 



have B 2 = B(l — t), and 

(j, = Lap + 2B(1 - t)p + 3B 3 p 2 . 

The critical point occurs when the first and second deriva- 
tives of the chemical potential are zero, giving two equa- 
tions for the critical value of t, t cp , and the critical density 

J- + 2B(l-t c P) + 6B 3 p cp = 



The second equation gives the density at the critical point 
straightaway, p cp = I/V6S3: the critical density does not 
depend on the value of the second virial coefficient but on 
that of the third virial coefficient. For our mixtures of ideal 
polymers and spheres, as the polymers are ideal their blobs 
do not repel each other, they only interact with each other 
via the spheres. Thus the third virial coefficient of the 
polymer blobs (at constant sphere chemical potential) is 
very small and the density at the critical point of demixing 
correspondingly high. 
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